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Abstract 

We investigate Hamiltonians with attractive interactions between pairs of 
fermions coupled to angular momentum J. We show that pairs with spin J 
are reasonable building blocks for the low- lying states. For systems with only 
a J = J m ax pairing interaction, eigenvalues are found to be approximately 
integers for a large array of states, in particular for those with total angular 
momenta I < 2j. For 1 = eigenstates of four fermions in a single- j shell we 
show that there is only one non-zero eigenvalue. We address these observations 
using the nucleon pair approximation of the shell model and relate our results 
with a number of currently interesting problems. 
PACS number: 21.60.Ev, 21.60.Fw, 24.60.Lz, 05.45.-a 
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I. INTRODUCTION 



Since pairing has proven to be important in atomic, nuclear, and condensed matter 
physics, pair truncation approximations to many body wavefunctions have been extensively 
studied. The first example is the seniority scheme, introduced by Racah [1,2], for the clas- 
sification of states in atomic spectra and later applied extensively in nuclear physics, where 
S pairs with angular momentum zero are related with a strong and attractive monopole 
pairing interaction. The second example is the interacting boson model (IBM) introduced 
by Arima and Iachello [3], where the low-lying excitations of complex even-even nuclei are 
described successfully by s bosons which correspond to correlated S nucleon pairs with 
angular momentum zero and d bosons which correspond to D nucleon pairs with angular 
momentum two. Again, the success of the IBM in nuclear physics partly comes from the 
validity of the pairing plus quadrupole-quadrupole force for effective interactions between 
valence nucleons. Monopole and quadrupole pairing are important as well in low and high 
temperature super conduct iviy in materials [4,5]. 

In this paper we investigate the general pair truncation approximation for fermions in a 
single- j shell. The examples explored may provide a clue as how to classify the states which 
come from the diagonalization of an attractive pairing interaction for which two fermions 
are coupled to an angular momentum J, 



Hj = - E A u A 
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M=-J 
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where [ ] J M means coupled to angular momentum J and projection M. Most of examples 
pursued in this paper have n — 4, where n is the nucleon number. 

We first point out in this paper that the low-lying eigenstates of (1) can be approximated 
by wavefunctions with pairs with angular momentum J only. We shall next show that a 



large array of eigenvalues of four nucleons in a single- j shell are asymptotic integers when 
J ~ Jmax = 2j — 1 and the total angular momentum, /, is not very close to I max = 4j — 6, 
and that this phenomenon originates from validity of the pair truncation scheme and special 
features of coupling coefficients. We shall finally prove that the pair Hamiltonian (1) has 
exactly one and only one non-zero eigenvalue for four fermion eigenstates with total angular 
momentum zero, 1 = 0. This sheds light on the problem of angular momentum zero ground 
state dominance in many-body systems interacting by random interactions [6]. 

II. COMPARISON OF THE PAIR APPROXIMATION TO THE SHELL MODEL 

Fig. 1 compares the exact ground state angular momentum / for four nucleons in a single- 
j shell interacting by the attractive pair Hamiltonian Hj with the angular momentum / of 
the ground state in the truncated space of pairs coupled to angular momenutm J only. 
We have examined all the cases up to J = 20 and j < 31/2 but here we show only two 
typical examples with J = 6 and 14. In the case of J = 0, the seniority scheme (S-pair 
approximation) produces the exact ground state. In case of J = 2 and n = 4, a /J-pair 
approximation is found to be always very good. When J > 2, the J-pair approximation of 
low- lying states is not perfect but always very reasonable. In Fig. 1, most of ground state 
angular momenta are correctly reproduced by the J-pair truncation. In Fig. lb), there are 
two exceptions in which the ground state is not correctly given by two J = 14 pairs: 2j = 19 
and 25. 

Even in those cases where the ground state angular momenta are not correctly predicted 
by the J-pair truncation, the low-lying state energies are reasonably reproduced (including 
the binding energies). As a "bad" example in which the ground state angular momentum 
is not reproduced, we show in Fig. 2 the case of four nucleons in a j = 25/2 shell with 
J = 14. The calculated levels using two J = 14 pairs are shown in the first column. The 
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next two columns are the shell model states obtained by diagonalizing the Hamiltonian in 
the full space. The states in the second column are the shell model states correponding to 
the pair truncation states in the first column. All the levels below + are included. One 
sees that the lowest four states 2+, 6+, 12+, and 10+ are well approximated by two J = 14 
pairs, although the ground state angular momentum is not correctly given by this J = 14 
pair approximation. The 0+ state is always precisely reproduced (as we shall see there is 
only one non-zero eigenvalue for / = states and that eigenstate is constructed by J pairs). 
The states with odd / are always outside the pair truncation space but their energies are 
quite high in all cases that we have checked. In Fig. 2, the 9 + is the lowest state with odd 
value of /. Below the 9+ there are 10 states with even values of / and most of which are 
nicely described by the J = 14 pair approximation. The angular momentum of those states 
for which excited energies are not well reproduced by two J pairs are labled by italic font 
in Fig. 2. 

For J = J max = 2j - 1, the I = I max = Aj - 6 or / = I max - 2 states are always 
the lowest. These two states may be constructed by pairs with angular momentum either 
Jmax or J max — 2. However, pairs with angular momentum J ma x — 2 do not present a good 
classification for other I states while those with angular momentum J max do. 

For n=3, the J-pair truncation describes the low-lying states precisely. For n=5 and 6 
cases that we have examined, the low-lying states are reasonably approximated by the J-pair 
truncation. We note without details that bosons with spin / exhibit a similar situation. It 
would be interesting to know the situation in more complicated systems. 

III. INTEGER EIGENVALUES 

We next report a very interesting regularity in the spectrum of the Hamiltonian (1) 
with J = J max = 2j — 1. The eigenvalues of most states with I < 2j — 3 are found to 
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be very close to integers corresponding to the number of pairs with angular momentum 
J max except for a very few eigenvalues. Taking four nucleons in a j = 31/2 shell as an 

example, the diagonalization of Hj (J = 30) gives "integer" eigenvalues for low I states 

all eigenvalues are 0, —1, and —2 to within a precision of 0.01 for all states with / < 22. For 
states with 22 < I < 52, these three "integer" eigenvalues continue to be valid except that 
7 eigenvalues which are not close to 0, —1 or —2 come in. These "non-integer" eigenvalues 
are very stable in magnitude for states with 22 < I < 52. The states with / > 53 are one 
dimensional, so the corresponding eigenvalues (saturate quickly with j) may be analytically 
derived [7]. These "integer" eigenvalues are best seen in case of J = J max and becomes less 
dominant for smaller J and the same single- j shell. 

To understand the validity of the pair approximation and the occurrence of these "inte- 
ger" eigenvalues we consider the pair basis of four nucleons 

1 /., + 



|rW2]/,M> 



V JV JiJ2;JiJ 2 



(A*t x A*t)£| ), 



(2) 



where Nj I ^j 2 . Ji j 2 is the diagonal matrtix element of the normalization matrix 

j j Jl 



(3) 



j j Ji 
J[ J> 2 I. 

In general this basis is overcomplete and the normalization matrix may have zero eigenvalues 
for a given /. 

The matrix elements of Hj are [8] 



(j i [J[j! 2 ]I,M\Hj\j 4 [J l J 2 }I,M) = 
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where Ji is a short hand notation of v / 2J 1 ~+T. There are two terms in 
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the 



second term is a nine-j coefficient which is usually much less than unity in magnitude, in 



particular when J is large and I not large (refer to Appendix A). Neglecting this nine- 
j symbol, the allowed states are \j 4 [J 1 J 2 ]I, M) « (A Jlt x ^ J2t )^ |0) for J 2 < J ± , and 
\j 4 [Ji Ji]I, M) ~ ^ ^A Jl ^ x A JlJ< yj |0), / even only, and the Hamiltonian matrix becomes 

4 [J X J 2 ]/, M\Hj\ 3 %h J 2 ]7, M) « 1 2 - == + <W)- ( 5 ) 

^i + (-iy5 Ju j 2 yfi + 

First of all we see that the matrix is diagonal which validates the pair approximation. Sec- 
ondly the eigenvalues are either 0, —1, or —2 with their corresponding wavefunctions having 
0, 1 or 2 pairs with angular momentum J, respectively. Therefore, the integer eigenvalues 
of many I states originate from both the special properties of these nine-j symbols and the 
validity of J-pair truncation. 

From the J-pair coupling scheme, the number of |j 4 [Ji J 2 ]I, M) with J 1 = J 2 = 2j — 1 
is 1+( ~ )J , and the number of \j 4 [JiJ 2 }I, M) with J x = 2j - 1, J 2 < J x and I < 2j - 1 is 
[ | ] (the largest integer not exceeding |). According to the above discussion, the number 
of states with eigenvalues ~ —2 is 1+ ^- > and the number of states with eigenvalues ~ — 1 
is [ | ]. This is confirmed in all cases with I < 2 j — 8. Eigenvalues not close to integers 
arise in states with 2j — 8 < I < 4j — 12. These "non-integer" eigenvalues are found to be 
almost the same for the 2j — 8 < I < Aj — 12 states; an understanding of this regularity is 
in progress. 

From a more general expression of Eq. (4), say, Eq. (5.8) of [9], we expect that the 
"integer" eigenvalues appear not only in even systems but also appear in odd-A systems. 
According to our numerical results, the pattern of "integer" eigenvalues also appear in the 
states with small I for j > 11/2 and n — 3, and for j > 23/2 and n = 5 etc. For cases 
with n = 3, a similar proof is readily obtained in terms of 6-j symbols. An explicit proof 
for more nucleons will be quite complicated. 
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IV. EXACT RESULTS FOR ANGULAR MOMENTUM ZERO STATES 



We now come to the last point of this paper by pointing out that there is only one 
non-zero eigenvalue for I = and n — 4. We define a new basis for the 1 = states for the 
pair Hamiltonian Hj 

N (o) 

Ifj,) = = M = 0) - b -4 [JJ]7 = M = 0), (J, + J), (6) 



JJ;JJ JiJi',JiJi 

\fj) = \f[JJ]I = M = 0). (7) 

The unnormalized states |j 4 Ji), Ji 7^ J, are orthogonal with respect to on/?/ |j 4 [JJ]J = M = 
0). From (4) one has 

N (0) N (0) 

(jV'A% ^\Hj\jVM^ o) = - ; JlJi;JJ JlJi;JJ , (8) 

r(0) at-(O) 



iV J(J(;JJ JV JiJi;JJ 

where J{, J± = 0, 2, • • • , 2j — 1. Using this formula, one easily confirms that all matrix 
elements of the Hamiltonian in the basis (6), (j 4 J{|iJj|j 4 Ji), J[ 7^ J, are zero. Therefore 
all the eigenvalues for n = 4 and 7 = are zero except for the single state with both pairs 
having angular momentum J and its eigenvalue is E'^ = —Nfj.jj. 

Since Hj is a negative definite operator, its eigenvalues will be negative or zero. From 
above we see that, for 1 = 0, there is only one state with a non-zero eigenvalue, 
Therefore one expects this eigenvalue to be the lowest in the spectrum because the eigen- 
values of / 7^ states are more or less scattered in many states, generally speaking. Thus 
the probability that the / = + is the lowest state of four fermions in a single-j shell in the 
presence of random two-body interactions is expected to be larger than the probability for 
all other angular momentum /, according to the empirical rule of Ref. [7]. For j < 31/2 
there are only two exceptions , j = 7/2 and 13/2. 

The sum rule of diagonal matrix elements [10] gives J2j — ~\ n { n ~ 1)Dq^ = 

—6D (j), where is the number of I = states and here n — 4. For n = 4 the number 



of states is D ( J) = [(2j + 3)/6] [11] which gives 1, 1, 1, 2, 2, 2, 3, 3, 3, • • • for 2j=3, 5, 7, 9, 
11, 13, 15, • • •, etc., regularly. Thus the staggering in the number of states is expected to 
be reflected in the staggering of the energy which was pointed out in [7,12] but without an 
explanation. 

V. SUMMARY 

In this paper we have shown that an attractive J-th pairing interactions favors pairs 
with angular momentum J in low-lying states of fermions in a single- j shell. Therefore, one 
may use pairs with angular momentum J as building blocks of wavefunctions of low-lying 
states. This is in contrast to repulsive pair interactions, used, for example, in the fractional 
quantum Hall effect [13-15], for which the pair truncation approximation is not valid. 

In addtion we discovered that the eigenvalues of states with low angular momentum / for 
pair Hamiltonians with J larger than j are approximately integers. We explain the origin of 
this fact for four nucleons in a single-j shell from the validity of pair truncation and special 
properties of nine-j symbols. We point out without details that the same holds for n = 3 
and j > 11/2 and for n = 5 (j < 23/2). 

Finally we pointed out that there is only one non-zero eigenvalue for / = and n = 4. 
This result, together with the empirical rule of ref. [7], provides a simple argument for the 
large probability of angular momentum zero states to be the lowest in energy for n = 4. 
Also for the first time, we demonstrated that the staggering of this probability for n = 4 vs. 
j is directly related to the staggering of number of / = states vs. j. 
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VII. APPENDIX A: SOME NINE-J SYMBOL FORMULAE 



These formulae are obtained in the following two steps: First, rewrite the nine-j in terms 
of six-j symbols, i.e., 



E(-f(2t + l) 



3 3 ri j j r 2 



n I t 



3 t s 2 




3 3 r 1 
3 3 r 2 
Si s 2 I 

and second, make use of the analytical formulas of six-j. Through these examples (though 
we are unable to get a universal formulae), one sees that the nine-j symbols in Eq. (3) are 
much less than unity and may be neglected in (4) when I is not very large. 



3 3 2 j - 1 

3 3 2 j - 1 

2 j - 1 2j - 1 

3 3 2 j - 1 

3 3 2 j - 1 

2 j - 1 2 j - 1 2 

3 3 2 j - 1 

3 3 2 j - 1 

2 j - 1 2j - 1 4 

J J 2 j - 1 

J J 2 j - 1 

2j - 1 2 j - 1 6 

J J 2 j - 1 

J J 2 j - 3 

2 j - 1 2 j - 3 2 



j(4j-3) [(2j-l)!] 2 
(4j-1)(4j-1)! ' 



j(8j 2 -6j-3) [(2j-l)!f 
(4j-1)(4j-3)(4j-1)! 



3j(2j + l)(4j 2 -3j-5) [(2j-l)!] 2 
(4j-1)(4j-3)(4j-5)(4j-1)! ' 



5j(j + l)(2j + l)(8j 2 - 6j - 21) [(2j - 
(4j - l)(4j - 3)(4j - 5)(4j - 7)(4j - 1)! 



> = 



3 (36 + j(4j - 9) (19 + 2j(~9 + 4j))) [(2j - 1)!] 
3(4j-3)(4j-5)(4j-1)! 



10 



3 J 2j - 1 
3 J 2j - 3 
2 j - 1 2 j - 3 3 



> = 



j(4j-3)(1 + 2j(4j-9))[(2j-1)!]- 
6(4j-5)(4j-l)! 



From the above formulae one sees that these nine-j symbols are proportional to 



[(2j-l)!] 



(4j-l)! ' 

and are very close to zero when j becomes considerably large. For example, the absolute 
values of these nine-j symbols are less than ~ 1CT 15 when j = 31/2. 



VIII. APPENDIX B: A NEW SUM RULE FOR A SIX-J SYMBOL 



From Eq. (3) one obtains 



Ar(° J ) ;JJ = 2 + 4(2J+l) 



j j J 
j j J 



> . 



Since -N { f) .jj is also the unique eigenvalue of the 1 = eigenstate of Hj, one has a sum 
rule that 



where n — 4, Dq 



(j) 



2j+3 
6 



even J 



[11]. One finally obtains that 



E (2-7 + 1) { 

even J 



j j J 

3 3 J 



> = < 



I if 2j = 3/c, 
0, if 2j + 2 = 3k, 
-\, if 2j -2 = 3k. 



(9) 
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FIGURES 

FIG. 1. Ground state angular momenta I for four fermions in a single- j shell for the the pair 
Hamiltonian Hj for J=6 in (a) and 14 in (b) as a function of j. The solid squares are ground 
state angular momenta obtained by truncating the space of states to those with two pairs with 
angular momentum J only, and the open circles are ground state angular momenta calculated by 
the diagonalizing the pair Hamiltonian in the full shell model space. 

FIG. 2. A comparison of low- lying spectra calculated with wavefunctions with two pairs with 
angular momentum J = 14 (the column on the left hand side) and by a diagonalization of the 
full space (the column in the middle and the column on the right hand side) for the case of four 
nucleons in a single- j (j = 25/2) shell. The middle column plots the shell model states which 
are well reproduced by the two J = 14 pairs, and the right column plots the shell model states 
which are not well reproduced by two J = 14 pairs. All the levels below 0^ in the full shell model 
space are included. One sees that the low-lying states with / = 2+, 6f, 12± , and 10^ are well 
reproduced. It is noted that the 0^ coupled by two J = 14 pairs is equivalent to that obtained by 
a full shell model diagonalization; refer to the text. 
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